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OVERVIEW 

Maintenance  or  restoration  of  physical 
aquatic  habitat  in  streams  during  critical 
periods  can  often  be  accommodated  with 
the  development  of  low-flow  channels, 
designed  to  concentrate  flows  and  increase 
channel  velocity  and  depth  during  low-flow 
periods.  A  procedure  for  the  first-order 
approximation  of  a  stable  low-flow  channel 
form  based  upon  physical  reasoning, 
empirical  evidence,  and  constraints  common 
to  low-flow  channel  projects  is  presented  in 
this  technical  note.  Topics  discussed 
include  simplifying  assumptions,  limitations, 
and  applicability. 

PLANNING 

Unaltered  natural  streams  support  healthy 
aquatic  ecosystems  due  largely  to  their 
habitat  complexity.  Alternate  pools  and 
riffles,  irregular  planforms,  eroding  banks, 
overhanging  vegetation,  boulders,  logs,  and 
variable  substrate  all  contribute  to  diversity 
of  habitat  within  and  along  the  channel. 
Natural  channels  are  scaled  to  provide  this 
diversity  for  a  wide  range  of  flow  conditions. 

In  contrast,  channels  modified  for  flood 
control  often  have  relatively  uniform  depths 
and  velocities  and  usually  contain  few 
instream  and  streambank  features  that 
provide  habitat.  The  same  situation  arises 
for  urban  streams  that  have  enlarged  due  to 
increased  runoff;  they  are  “overfit”  during 
baseflow  conditions  and  display  small, 
featureless  channels.  This  lack  of  diversity 
in  modified  and  enlarged  channels 
negatively  influences  the  aquatic  ecosystem. 


Attaining  habitat  diversity  in  such  channels 
requires  the  incorporation  of  instream  and 
streambank  features  and  cover  devices  that 
create  depth,  velocity,  and  cover  diversity. 
These  features  are  most  effective  when  they 
contribute  to  natural  and  diverse  channel 
conditions  during  low  flow  but  have  little 
impact  upon  the  ability  of  the  channel  to 
convey  flood  flows  and  sediment. 
Furthermore,  these  features  must  possess 
the  necessary  stability  to  withstand  the  high 
velocities  typically  associated  with  modified 
channels.  A  successful  design  concept  that 
improves  habitat  conditions  but  has 
negligible  effects  upon  flooding  is  the 
combination  of  streambank  and  instream 
features  to  form  a  diverse,  stable,  low-flow 
channel  (Figure  1). 


Figure  1.  Low-flow  channels  can 
maintain  natural  stream  characteristics 
within  an  enlarged  floodway 
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Restoration  of  physical  aquatic  habitat  in 
streams  during  critical  periods  can  often  be 
accomplished  with  low-flow  channels 
designed  to  concentrate  flows,  increase 
channel  velocity  and  depth,  and  generally 
maintain  the  characteristics  of  a  natural 
channel  within  an  enlarged  flood  control 
channel.  Structural  measures  such  as  bank 
protection,  flow  deflectors,  and  sills  are 
generally  required  to  stabilize  the  channel 
and  provide  habitat  diversity.  In  addition  to 
the  usual  design  requirements  for  these 
structures,  a  structural  layout  is  needed  that 
allows  the  stream  to  adopt  a  form  consistent 
with  the  project  objectives.  Designers  must 
translate  physical  habitat  objectives  into 
hydraulic  and  geomorphic  criteria,  which,  in 
turn  are  used  to  select  structural  measures 
(Fischenich  et  al.  1994).  Unfortunately,  no 
procedure  has  been  proposed  to  accomplish 
this. 

In  this  technical  note,  a  stable  channel  is 
defined  as  a  channel  that,  over  an 
engineering  time  scale,  retains  its  planform, 
cross-sectional  geometry,  and  grade  with 
allowance  for  temporal  variability  about  a 
mean  that  displays  no  trends.  Existing 
techniques  for  predicting  the  hydraulic 
geometry  of  stable  alluvial  channels  can  be 
placed  into  two  categories:  1)  empirical 
regime  equations  derived  from  regression 
analyses  of  observed  channel  geometry, 
and  2)  analytical  models  that  attempt  to 
model  rivers  on  a  rational  basis  using 
theoretical  considerations.  Knighton  (1984), 
Chang  (1988),  and  Mueller  and  Dardeau 
(1990),  among  others,  present  summaries  of 
many  of  these  techniques.  The  early  regime 
equations  generally  relate  some  two- 
variable  combination  of  velocity,  width, 
depth,  slope,  area,  hydraulic  radius, 
perimeter,  or  discharge  as  a  power  function 
and  can  be  generalized  in  the  form  of 
Equations  1  through  3,  where  by  continuity, 
ack=1  and  b+e+m=1 :2 

W  =  aQb  (1) 

D  =  cQe  (2) 
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V  =  kQm  (3) 

Analytical  procedures  include  the  method  of 
maximum  permissible  velocity,  the  critical 
shear  stress  method,  and  the  so-called 
“rational”  regime  equations  that  use  tractive 
force  theory,  resistance  and  sediment 
transport  laws,  minimum  stream  power, 
maximum  sediment  transport  capacity, 
similitude,  bank  stability  analyses,  and  other 
concepts  to  overcome  the  fact  that  alluvial 
channels  have  more  degrees  of  freedom 
than  available  relations  for  solution. 
Wargadalam  (1993)  showed  that  the 
existing  body  of  literature  on  the  subject 
could  be  reduced  to  no  fewer  than  176 
equations  describing  depth,  width,  velocity, 
and  slope  in  the  form: 

W,D,V,orS0=aiQb'dsc'Qse'  (4) 

Despite  the  considerable  effort  expended  in 
formulating  these  many  approaches  to 
stable  channel  geometry  prediction,  none 
are  suitable  for  the  development  of  stable 
low-flow  channel  designs.  The  reasons  for 
this  include:  1)  no  consideration  of  habitat 
requirements,  2)  empiricism,  3)  applicability 
only  to  straight  channels,  and  4)  theoretical 
flaws.  Perhaps  the  greatest  limitation  of 
many  of  these  techniques  is  that  they 
consider  width,  depth,  and  slope  to  be  the 
only  degrees  of  freedom  and  do  not  account 
for  the  full  dynamic  nature  and  multi¬ 
dimensionality  of  fluvial  systems.  A  minimal 
description  of  channel  morphology  requires 
that  its  grade,  planform,  cross-sectional 
geometry,  and  some  measure  of  channel 
resistance  be  parametrically  defined. 
Channel  grade  can  be  described  by  the  bed 
slope  S0  and  resistance  can  be  described  by 
the  friction  slope  Sf.  Channel  planform 
requires  at  least  two  parameters,  e.g., 
sinuosity  and  meander  arc  length.  To 
uniquely  define  the  cross-sectional  geometry 
of  a  channel  requires  that  at  least  one 
variable  in  addition  to  the  width  and  depth 
be  specified.  Examples  of  a  third  variable 
are  maximum  depth,  bank  slope,  and  a 
channel  shape  factor.  Figures  2  and  3  show 
the  parameters  used  in  describing  channel 
geometry  and  planform  in  this  paper. 
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See  page  9  for  mathematical  notation. 
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L  meander  wavelength 
M f  meander  arc  length 
w  average  width  at  bankful l  discharge 
Ma  meander  amplitude 
rc  radius  of  curvature 
0  arc  angle 

Figure  2.  Channel  planform  definition 
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Figure  3.  Channel  cross  section  definition 
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PROCEDURE  FORMULATION 

Fundamental  Considerations 
Both  Cartesian  and  curvilinear  coordinate 
systems  are  used  in  the  formulation  of  the 
procedure.  The  principal  axis  x  in  the 
Cartesian  system  defines  the  centerline  of 
the  meandering  pattern  in  the  downstream 
direction  for  planform  analysis;  it  is  aligned 
with  the  stream  and  is  horizontal  for  cross 
section  and  grade  analyses.  The  y  and  z 
axes  alignments  are  with  the  vertical  and 
horizontal,  both  normal  to  the  flow.  In  the 
curvilinear  system,  the  sinuous  axis  s 
follows  the  centerline  of  the  stream,  and  the 
transversal  direction  remains  orthogonal  to 
the  principal  axis  in  the  horizontal. 

Analytical  determination  of  the  hydraulic 
geometry  for  alluvial  channels  can  be 
accomplished  only  if  applicable  physical  and 
empirical  relations  are  sufficient  to  describe 
the  unknowns  or  degrees  of  freedom.  Two 
characteristics  of  low-flow  channels  help  to 
satisfy  this  requirement.  First,  habitat 
analyses,  based  on  biological  criteria,  are 
used  to  establish  a  minimum  flow  depth  and 
an  optimum  velocity  for  the  channel. 

Second,  instream  and  streambank 
structures  are  typically  used  to  stabilize  the 
channel,  provide  cover  and  substrate,  and 
generate  diversity  in  velocity.  Given  these 
constraints,  physical  reasoning,  geometric 
analyses,  and  a  few  semiempirical  relations, 
a  solution  is  possible. 

A  fundamental  assumption  in  the  application 
of  this  procedure  is  that  sand  bed  channels 
tend  toward  a  meandering  planform  that  can 
be  described  by  a  sine-generated  curve. 
Langbein  and  Leopold  (1966)  used  the 
theory  of  minimum  variance  to  show  that  the 
most  probable  path  of  a  meandering  stream 
was  a  sine  curve.  Julien  (1985)  showed  that 
the  fundamental  shape  of  meandering 
planforms  could  be  evaluated  based  on  the 
separation  of  the  boundary  layer  near  the 
inner  bank  and  the  rate  of  energy 
dissipation.  Julien  (1985)  demonstrated  that 
among  all  possible  geometries,  with  the 
exception  of  a  straight  channel,  boundary 


separation  and  the  rate  of  energy  dissipation 
are  minimized  by  a  simple  sine-generated 
function.  Approaches  based  upon  the 
Kinoshita  equation,  Von  Schelling’s  curve, 
and  Fargue’s  spiral  can  also  be  reduced  to  a 
form  of  sine-generated  curve.  The  sine¬ 
generated  function  has  thus  been  shown  to 
fit  the  shape  of  meanders  on  stable  rivers 
quite  well  and  is  the  most  convenient 
expression  to  use.  The  corresponding 
meandering  pattern,  with  reference  to 
Figure  2,  can  be  written  as: 

dQ  .  2 ns 

—  =  -a,  sin -  (5) 


A  second  fundamental  assumption  is  that 
the  flow  is  approximately  uniform  in  the 
streamwise  direction.  This  assumption 
allows  the  bed,  water,  and  friction  slopes  to 
be  equated  and  permits  the  use  of  a  uniform 
rating  equation.  The  final  assumption 
relates  to  sediment  continuity,  which  must 
be  preserved  to  meet  the  stable  channel 
definition.  The  proper  approach  to  ensure 
that  continuity  is  met  is  to  perform  both 
longitudinal  and  transversal  sediment 
transport  analyses.  A  simplified  approach  is 
used,  wherein  continuity  is  approximated  by 
setting  the  bed  shear  stress  at  the  critical 
value  for  the  mean  bed  sediment  size. 

Procedure 

Because  low-flow  channels  are  developed 
for  environmental  purposes,  the  first  step  is 
to  conduct  hydrologic  analyses  to  establish 
the  target  or  design  discharge.  Instream 
Flow  Incremental  Methodology  (IFIM) 
analyses,  habitat  suitability  curves,  or  other 
biological  criteria  are  then  used  to  establish 
target  values  of  velocity  and  depth. 

Realizing  optimum  values  for  both  the  depth 
and  velocity  objectives  is  unlikely;  therefore, 
one  must  be  established  as  a  critical 
parameter.  The  most  common  case,  in 
which  a  critical  minimum  depth  is 
established  and  an  optimum  velocity  target 
identified,  is  used  here.  Simple 
modifications  to  this  procedure  will 
accommodate  other  cases,  such  as  a 
maximum  velocity  and  an  optimum  depth. 
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Grade:  The  criterion  of  beginning  of  motion 
establishes  a  threshold  corresponding  to  a 
condition  of  maximum  allowable  shear.  The 
ratio  of  forces  on  an  individual  sediment 
grain  on  a  sloping  bed  is  evaluated  to  define 
the  dimensionless  shear  stress  called  the 
Shields  number: 


_  yRhS / 

(y,-yK 


(6) 


The  critical  value  of  the  Shields  number  t*c 
corresponding  to  the  beginning  of  motion, 
was  determined  by  Shields  (1936).  For 
particle  Reynolds  numbers  greater  than  200, 
T*c  =  0.05.  The  widely  accepted  Shields 
diagram  is  used  as  a  basis  for  many  stable 
channel  design  techniques  despite 
complexities  encountered  in  its  application 
to  irregular  channel  geometries.  Given  the 
depth  constraint  and  sediment  size 
measurements,  the  maximum  allowable 
friction  slope  for  longitudinal  stability  can  be 
computed  from  Equation  6  by  setting  the 
hydraulic  radius  Rh  equal  to  the  minimum 
depth  and  setting  t*c  =  0.05.  For  sands  and 
gravels  of  specific  gravity  G  =  2.65,  this 
relation  reduces  to: 


Sf=  0.0825^  (7) 

For  uniform  flow,  channel  grade  can  be 
equated  to  the  friction  slope,  and  either  can 
be  related  to  the  sinuosity  of  the  low-flow 
channel  and  measurement  of  the  bed  slope 
of  the  parent  channel  Sc: 


which  can  be  incorporated  into  the  analysis 
as  discussed  later. 

Planform:  With  the  allowable  friction  slope 
computed  by  Equation  7  and  a  measured 
thalweg  slope  of  the  parent  channel, 
Equation  8  is  used  to  establish  the  ratio  of 
the  meander  wavelength  to  the  arc  length. 
Again,  assuming  a  sine-generated  curve  for 
the  planform  shape  and  a  constant  velocity 
along  the  meander  length,  Equation  5  can 
be  integrated  for  0: 


^  Ma7  27 is 

0  = - cos - 1-  c 

2nV  M 


(9) 


in  which  the  maximum  angle  of  the 
streamline  with  the  downstream  valley  0m 
corresponds  to: 


0 


m 


Ma2 

2jiV 


(10) 


The  meander  wavelength  is  computed  from 
the  following  relation: 


M 

k  =  jcos0rA  (11) 

o 

Defining  a  nondimensional  curvilinear 
distance,  s’  =  s/M,  and  rearranging  with 
Equations  9  and  10  gives: 

*1  1 

—  =  |  cos[0  m  cos(27u')]  ds'  (12) 
M  0 


While  this  must  be  the  reach-wise  average 
bed  slope  of  the  low-flow  channel,  variation 
about  this  average  to  accommodate  pool- 
riffle  sequences  is  possible  if  the  resultant 
reaches  have  the  requisite  transport 
capacities  to  ensure  no  aggradation  or 
deposition.  This  first  approximation  of  the 
channel  grade  ignores  transversal  stability, 


The  sinuosity  is  defined  as  the  ratio  M/k  and 
increases  gradually  with  0m.  Equation  12 
has  been  integrated  numerically  and  the 
results  are  shown  in  Figure  4.  Figure  4  (or 
Equation  12)  can  be  used  to  determine  the 
maximum  angle  of  the  streamline  0m.  For  a 
given  friction  slope  and  the  assumption  of  a 
sine-generated  curve,  the  shape,  sinuosity, 
and  maximum  angle  of  the  streamline  are 
defined.  The  scale  must  still  be  determined 
even  though  the  shape  of  the  planform  is 
known. 
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The  radius  of  curvature  varies  along  the 
bend  as  a  cosecant  function,  and  the 
minimum  radius  of  curvature  rm  is  reached 
when  the  cosecant  function  equals  unity. 
The  ratio  of  the  wavelength  to  the  minimum 
radius  of  curvature  X/rcm  is 


X  _  2%QmX 

r~  M 

m 


(13) 


The  meander  amplitude  Am  as  defined  in 
Figure  2  is  evaluated  analytically  by  the 
following  integral: 


M/4 

Am  =  2  j  sin0  ds  (14) 

o 

The  ratio  of  wavelength  X  to  the  meander 
amplitude  Am  is 


cos[0  m  cos(27u' 

0 


2  j  sin[0  m  cos(27U')]  ds' 

o 

Equations  13  and  15  have  been  integrated 
numerically  and  are  also  shown  in  Figure  4. 
The  ratio  X/Am  is  a  function  of  0m  that 
increases  rapidly  when  0m  exceeds  90°  and 
reaches  a  value  of  3.25  at  the  meander 
cutoff  (0m  s  125°).  The  relations  expressed 
in  Equations  13  and  15  further  define  the 
planform  shape,  but  an  additional  relation  is 
needed  for  the  size.  One  relation  can  be 
formulated  on  the  basis  of  a  limiting  radius 
of  curvature  to  ensure  cross-sectional 
stability.  This  procedure  is  tedious  and 
requires  an  analysis  of  the  transversal  shear 
stress  distribution  and  sediment  transport 
along  the  meander  path.  In  the  case  of  low- 
flow  channel  design,  a  simpler  relation  is 
available  in  that  the  meander  amplitude  is 
constrained  by  the  width  of  the  parent 
channel.  The  latter  approach  has  been 
adopted. 


0  20  40  60  80  100  120 


A  % 

Figure  4.  Planform  relations 


6 


ERDC  TN-EMRRP-SR-1 9 


Low-flow  channels  are  typically  constructed 
by  placing  longitudinal  bank  protection,  flow 
deflectors,  and  sills  within  a  parent  channel. 
The  width  of  the  parent  channel  at  the  low- 
flow  channel  elevation  defines  a  maximum 
meander  amplitude  with  no  erosion  of  the 
parent  banks.  This  is  only  an  upper  limit, 
and  any  0  <  Am  <  Wc  will  meet  the  shape 
requirements.  From  a  practical  standpoint, 
however,  constructability  and  habitat 
benefits  are  optimized  with  the  larger 
amplitudes.  The  amplitude  should  be  set  at 
some  fraction  (perhaps  90  percent)  of  the 
parent  channel  width,  minus  the 
encroachment  of  longitudinal  revetment 
structures.  With  this  value  of  Am,  either 
Figure  4  or  Equations  12  and  15  can  be 
used  to  determine  the  meander  arc  length 
and  wavelength.  The  channel  planform  can 
then  be  laid  out,  and  structures  can  be 
selected  and  located  in  the  normal  fashion. 


Cross  Section:  Although  habitat  criteria 
were  used  to  establish  a  target  velocity,  a 
rating  equation  is  used  to  determine  the 
mean  velocity  for  the  defined  channel.  The 
Manning-Strickler  Equation  can  be  used  to 
find  mean  velocity  V: 


V  = 


21.1 

d w 

u50 


2/3  cl/2 
/ 


Dz,iS 


(16) 


This  mean  velocity  must  be  evaluated  on 
the  basis  of  habitat  considerations.  If  it  is 
unacceptable,  the  depth  must  be  adjusted, 
and  the  entire  procedure  must  be  iterated 
until  an  acceptable  value  of  velocity  is 
obtained.  When  the  velocity  is  within 
acceptable  limits,  the  principle  of  continuity 
is  applied  to  determine  the  cross-sectional 
area: 


A  = 


Q 

V 


(17) 


An  infinite  number  of  cross-sectional  shapes 
can  be  developed  within  the  area  and  depth 
constraints  given.  However,  naturally  stable 
sand  bed  channels  tend  to  have  certain 


shapes  with  respect  to  their  position, 
commonly  trapezoidal  in  the  crossings  and 
a  skewed  parabola  in  the  bendways. 
Investigators  have  described  these  shapes 
in  terms  of  transversal  sediment  transport, 
force  analyses,  and  critical  transversal 
shear  stress  distribution.  From  a  practical 
standpoint  these  complex  cross  sections 
would  be  difficult  to  construct,  and  the  width 
varies  little  from  that  of  a  trapezoidal 
section.  As  a  first  approximation,  a 
trapezoidal  section  with  three-to-one  side 
slopes  should  be  used  to  establish  the 
channel  bottom  width: 

W  =  —  -  3D  (18) 

D 


DESIGN  REFINEMENT  AND 
IMPLEMENTATION 

The  success  of  this  technique  on  sand-bed 
channels  is  virtually  nonexistent  due  to  the 
high  transport  rates  of  sediment  under 
nearly  all  flow  regimes.  Success  in  gravel- 
bed  streams  is  only  slightly  higher  than  for 
sand-bed  streams  (Fischenich  1993). 
Therefore,  accomplishment  of  project 
objectives  can  be  significantly  improved  by 
using  structural  measures  such  as  bank 
protection,  flow  deflectors,  and  sills  to 
stabilize  the  planform,  cross-sectional 
geometry,  and  bed  profile  of  the  channel.  In 
addition  to  the  usual  design  requirements 
for  these  structures,  the  designer  must 
select  an  appropriate  structural  layout  with 
the  objective  of  causing  the  stream  to  adopt 
a  form  that  is  consistent  with  the  project 
objectives.  Projects  are  too  often 
formulated  with  little  or  no  thought  given  to 
the  placement  of  structures;  as  a 
consequence,  both  structural  and  channel 
stability  are  compromised.  With  no 
assurance  of  the  future  form  of  the  channel, 
projected  benefits  cannot  be  assured. 
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APPLICABILITY  AND 
LIMITATIONS 

A  procedure  for  the  first-order 
approximation  of  a  stable  low-flow  channel 
form  has  been  developed  based  upon 
physical  reasoning,  empirical  evidence,  and 
constraints  common  to  low-flow  channel 
projects.  More  sophisticated  analyses 
would  be  warranted  to  design  the  channel 
as  nearly  like  its  stable  form  as  possible. 
However,  this  procedure  yields  a 
reasonable  estimate  of  the  planform  and 
grade  of  the  channel  and  maximizes 
probability  that  natural  transport  processes 
will  quickly  re-form  the  cross  section  without 
impacts  to  the  channel  stability. 

Many  of  the  assumptions  used  in  the 
formulation  of  this  procedure  are  not  always 
valid.  For  example,  the  Manning-Strickler 
Equation  is  assumed  to  apply,  which  further 
implies  that  the  sediment  is  in  the  sand  size 
range  and  that  the  flow  is  approximately 
uniform.  Only  the  longitudinal  shear  stress 
is  considered  in  this  simplified  analysis.  A 
detailed  analysis  should  also  consider  the 
transversal  stress  distribution,  which  usually 
results  in  a  decrease  in  the  allowable  friction 
slope.  Finally,  section  stability  under  a 
range  of  discharges  has  not  been 
established.  Such  a  determination  would 
require  that  friction  slope  decrease  with 
increasing  depth  and  that  transport  capacity 
through  the  section  match  the  inflowing 
load. 

Despite  these  limitations,  this  method  as 
presented  yields  a  reasonable  first-order 
approximation  of  the  stable  channel  form. 
This  “stable”  channel  configuration  can  be 
further  evaluated  using  more  sophisticated 
techniques.  The  proposed  procedure  yields 
a  close  approximation  of  natural 
occurrences  and  thus  allows  the  stream  to 
make  minor  adjustments  that  conform  to  its 
newly  introduced  geometry. 
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NOTATION 

A  cross-sectional  area  of  the  channel 

Am  meander  amplitude 

a, ai,...  empirical  coefficients 

b, bi,...  empirical  coefficients 

c, ci,...  empirical  coefficients 
ds  sediment  size 

D  average  flow  depth 
Dc  average  depth  in  parent  channel 
d50  representative  bed  sediment  size 
e,k,m  empirical  coefficients 
Q  total  discharge 
Qs  sediment  load 
rm  minimum  radius  of  curvature 
Rh  hydraulic  radius 
M  wavelength  path  distance 
s  curvilinear  longitudinal  distance 
s’  dimensionless  curvilinear  distance 

Sf  friction  slope 

S0  bed  slope 

Sc  bed  slope  of  parent  channel 
V  cross-sectional  average  downstream 

velocity 

W  channel  width 

Wc  parent  channel  width 

w  curvilinear  transversal  coordinate 

x  Cartesian  longitudinal  coordinate 

y  vertical  coordinate 

z  Cartesian  transversal  coordinate 

y  unit  weight  of  water 

ys  unit  weight  of  sediment 

X  meander  wavelength 

0  angle  of  streamline  with  x 

0m  maximum  angle  of  streamline  with  x 

x*  dimensionless  Shields  number 

x*c  critical  Shields  number 
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